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We investigate several models described by real scalar fields, searching for topological defects, 
and investigating their linear stability. We also find bosonic zero modes and examine the thermal 
corrections at the one-loop level. The classical investigations are of direct interest to high energy 
physics and to applications in condensed matter, in particular to spatially extended systems where 
fronts and interfaces separating different phase states may appear. The thermal investigations show 
that the finite temperature corrections that appear in a specific model induce a second order phase 
transition in the system, although the thermal effects do not suffice to fully restore the symmetry 
at high temperature. 
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I. INTRODUCTION 



Models described by real scalar fields in (1,1) space-time dimensions are among the simplest systems that support 
topological solutions. Usually, these topological solutions are named kinks, which are classical static solutions of the 
equations of motion, and the topological behavior is related to the asymptotic form of the field configurations, which 
has to differ in both the positive and negative space directions. To ensure that the classical solutions have finite 
energy, one requires that the asymptotic behavior of the solutions is identified with minima of the potential that 
defines the system under consideration, so in general the potential has to include at least two distinct minima in order 
for the system to support topological solutions. 

We can investigate real scalar fields in (3, 1) space-time dimensions, and now the topological solutions are named 
domain walls. These domain walls are bi-dimensional structures that carry surface tension, which is identified with the 
energy of the classical solutions that spring in (1, 1) space-time dimensions. The domain wall structures are supposed 
to play a role in applications to several different contexts, ranging from the low energy scale of condensed matter 
[l, 2, 3, 4] up to the high energy scale required in the physics of elementary particles, fields and cosmology ^, 6, 7]. 

There are at least two classes of models that support kinks or domain walls, and we further explore such models 
in the next Sec. In the first class of models one deals with a single real scalar field, and the topological solutions 
are structureless. Examples of this are the cfA and (ffi models 0- In the second class of models we deal with systems 
defined by two real scalar fields, and now one opens two new possibilities: domain walls that admit internal structure 
|El^llO>llBllJillMlIl|> and junctions of domain walls, which appear in models of two fields when t he p otential contains 
non-collinear minima, as recently investigated for instance in Refs.flllllllllilllllMlllllllHlMlllllllll. 

There are other motivations to investigate domain walls in models of field theory, one of them being related to 
the fact that the low e nerg y world volume dynamics of branes in string and M theory may be described by standard 
models in field theory |2M l29l y^l ■ Besides, one knows that field theory models of scalar fields may also be used to 
investigate properties of quasi-linear polymeric chains, as for instance in the applications of Refs. 'sil l3^ IssL l3^ . to 
describe solitary waves in ferroelectric crystals, the presence of twistons in polyethylene, and solitons in Langmuir 
films. 

Domain walls have been observed in several different scenarios in condensed matter, for instance in ferroelectric 
crystals 0], in one-dimensional nonlinear lattices |35|, and more recently in higher spatial dimensions - see |36j | and 
references therein. The potentials that appear in the models of field theory that we investigate in the present work 
are also of interest to map systems described by the Ginzburg-Landau equation, since they may be used to explore 
the presence of fronts and interfaces that directly contribute to pattern formation in reaction-diffusion and in other 
spatially extended, periodically forced systems [3. 13. l3^ l3^. "W. 'w\ . 

In the present work, in Sec. ^1 we review some known facts on kinks and domain walls, and there we also introduce 
other results. In particular, we present some new models described by one and by two real scalar fields. The 
investigations follow in Sec. IIIII where we search for the topological structures that generate kinks and walls. We 
reserve for Sec. II VI the study of stability of the solutions that we found in the former section. We investigate the finite 
temperature effects in Sec. where we examine the effective potential for a specific model, which engenders three 
minima at the classical level. We end the work in Sec. IVII where we comment on applications to nonlinear science 
and briefly review the results of the present investigation. 
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II. GENERAL CONSIDERATIONS 

In this work we are interested in field theory models that describe real scalar fields and support topological solutions 
of the Bogomornyi-Prasad-Sommerfield (BPS) type [4ll li^. In the case of a single real scalar field 0, we consider 
the Lagrange density 

/:=^9a</'9"0-V(0) (1) 

Here is the potential, which identifies the particular model under consideration. We write the potential in the 

form V{(j)) — (1/2) W^, where W = W{(j)) is a smooth function of the field (j>, and = dW/d(j). In a supersymmetric 
theory W is the superpotential, and this is the way we name W in this work. 
The equation of motion for (j) — (j){x, t) has the general form 



and for static solutions we get 



^ = W^W^^ (3) 



It was recently shown in Refs. |43l |44| that this equation of motion is equivalent to the first order equations 

if one is searching for solutions that obey the boundary conditions lim-c^-oo '/'(a;) = and \mi.x^-oa{d(j)/ dx) — 0, 
where (jii is one among the several vacua {^i, 02, ■•■} of the system. In this case the topological solutions are BPS (+) 
and anti-BPS (-) solutions. Their energies get minimized to the value = |AWij|, where I^Wij — Wi — Wj, with 
Wi standing for W{(t)i). The BPS and anti-BPS solutions are defined by two vacuum states belonging to the set of 
minima that identify the several topological sectors of the model. 

In the case of two real scalar fields and x the potential is written in terms of the superpotential, in a way such 
that V{<j), x) = (1/2) + (1/2) W^. The equations of motion for static fields are 



fx 

dx'^ 



W^W^^ + W^W^4, (5) 
W^W^^ + W^W^^ (6) 



which are solved by the first order equations 



dx dx 



Solutions to these first order equations are BPS (+) and anti-BPS (-) states. They solve the equations of motion, and 
have energy minimized to V'^ — \/SW^^\ as in the case of a single field; here, however, AVF*-' — Vl^(?!'j,Xi) ~^{4>3TXj)i 
since now we need a pair of numbers (0i,Xi) to represent each one of the vacuum states in the system of two fields. 
In the plane (0, x) we may have minima that are non coUinear, opening the possibility for junctions of defects. In 
the case of two real scalar fields, we can find a family of first order equations that are equivalent to the pair of second 
order equations of motion, but this requires that Wipip -f W^x = 0, in the case of harmonic superpotentials t43- 44]. 

We now turn attention to kinks and domain walls. Perhaps the most known example of this is given by the 4)^ 
model, defined by the potential !/((/)) = (1/2) {(jp - if. Here we are using natural units, and dimensionless fields 
and coordinates. In this model the domain wall can be represented by the solution 4>s{x) = ±tanh(x). The above 
potential can be written with the superpotential W{4)) — (p — 0'^/3, and the domain wall is of the BPS or anti-BPS 
type. The wall tension corresponding to the BPS wall is tg = 4/3. 
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We can also find structureless domain walls in other models, for instance in the 0^ model, which is described by 
the potential V^(0) — (1/2)0^ (0^ — 1)^. Here we have W^(0) = (1/2)0^ — (l/4)(/)'*, and the wall configurations are 
also of the BPS type, and are given by 0^ = (1/2)[1 ± tanh(a;)]. The wall tension is now tg — 1/4. This potential was 
investigated for instance in Ref. '45^1 . 

We can build another class of models, which is described by two real scalar fields. In this case the domain walls may 
engender internal structure. This line of investigation follows as in Refs. 0, 0, 0| and we illustrate such possibility 
with the system defined by the potential 

r(l + 2r)0\2 (9) 

where the parameter r 7^ is real. This model was first investigated in Ref. This potential follows from the 

superpotential x) ^ (1/3)'/'^ — fcj)x^, and the system supports the two-field solutions (f>{x) = tanh(2ra;) and 
= asech(2rx), with = 1/r — 2. These solutions are BPS solutions, and now the parameter r is restricted to the 
interval r G (0, 1/2). The limit r ^ 1/2 lead us to the one-field solution (f){x) = tanh(a::) and xi^) = 0- The two-field 
solutions obey (jp +x'^ /a^ — 1, which describes an elliptic arc connecting the two minima (±1,0) of the corresponding 
potential in the ((/>, x) plane. The one-field solutions represent standard domain walls, while the two-field solutions 
may be seen as domain walls having internal structure: the vector (</>, x) in configuration space describes an straight 
line segment for the one-field solution, and an elliptic arc for the two-field solution, resembling light in the linearly 
and elliptically polarized cases, respectively. The same solutions appear in condensed matter, in the anisotropic XY 
model used to describe ferromagnetic transition in magnetic systems, and there they are named Ising and Bloch walls, 
respectively - see for instance Ref. 0, chapter 7. 

In the present work we are especially interested in three other models. They are defined by potentials that describe 
a single real scalar field 

^1(0) = 101 + ^0' (10) 
^2(0) = + V (11) 
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and two real scalar fields 



1^3(0, X) = ^{1 - \4>\f - r{l - mx^ + 

lr\^ + 2r^<P\^ (12) 

The first model was investigated in Refs. [HO] to model an exactly soluble linear chain system. It has also been 
investigated more recently in Ref. jsj] and also in Ref. with different motivations. The other models are new, 
and some of their classical and quantum features will be examined below. 

We notice that the above models can be described by the following superpotentials 

WM) = (13) 

W2{cl>) = l^f" -\\4>\<f' (14) 

and also 

VF3(0,x) = </>-i|</.|(/.-r</)x' (15) 

The presence of superpotentials help simplifying the investigation, since the BPS solutions satisfy first order differential 
equations, which are simpler to solve compared to the equations of motion. 

Our interest in the first model is directly related to the fact that it is very much similar to the standard (j)^ model. 
Thus, we explore its classical solutions, to see how similar they are to that of the 0'* model. The same for the 
second model, which was invented to be similar to the 4>^ model, and also for the last model, which is an extension 
of the first model to the case of two fields, in a way similar to the model defined by Eq. ||5J), as done in Ref. [isj . 
These potentials are also of interest in condensed matter, in the investigation of spatially extended, periodically forced 
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systems, governed by the Ginzburg-Landau equation. Explicitly, the Ginzburg-Landau equation for a single real order 
parameter A — A{x, t) can be written in the form 

dA _ d^A dU 

'dt~'d^~dA ^^^^ 

where U — U{A) is a function of A, which may be identified with the potential of the field models that we are using 
in the present work. These models may then provide diverse scenarios for fronts and interfaces to travel according to 
the Ginzburg-Landau dynamics. 

III. TOPOLOGICAL SOLUTIONS 

In this section we search for BPS solutions in the models introduced in the former Sec. |nj We investigate the first 
order equations corresponding to each one of the three models separately. 

A. Model 1 

This model is defined by the superpotential of Eq. H13() . In this case there are two singular points, that represents 
the two minima of the potential in Eq. I|1U|) . These minima are wi = 1, U2 = — 1. The equation of motion for (j) = (f){x) 
is 

We are searching for topological solutions, so we can consider instead of the second order differential equation the two 
first order equations 

(IS) 



They support the BPS states 



After inspecting the explicit form of these BPS solutions, we notice that they have the standard kink-like profile, as 
in the (p'^ model, as we show in Fig. 1. Also, they arc linearly stable, minimizing the energy to ii = 1, below the value 
tg = 4/3 that appears in the t/)^ model. 

The interesting feature of this model is that at the classical level it is similar to the model, although the potential 
goes up to the second order power in the field. We illustrate the present model and its topological solutions in Fig. 1, 
where we also plot the topological state of the standard (/)^ model to offer a visual comparison between the BPS states 
of the two models. There we notice that the BPS states of model 1 are thicker than the corresponding states of the 
(j)'^ model. Also, we notice that the tension of these BPS states can be written as ti = {3/4)ts. 

We notice that this model may be used as an alternative to the 0'* model, in applications in condensed matter. 
For instance, it may represent the motion in the polyacetylene (PA) chain, describing single-double or double-single 
bound alternation in the PA chain, but this is out of the scope of the present work. 

B. Model 2 

Now we explore the second model, described by the potential of Eq. Hll|) . It can be written in terms of the 
superpotential of Eq. (|14|) . This model contains three minima, at the values ui = 1, W2 = — 1, W3 = 0. It supports two 
topological sectors of the BPS type. These sectors are degenerate, and the wall tension is given by ^2 = 1/6. The 
BPS solutions satisfy 



(20) 
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The solutions are 

(j)1i^{x) = ±i[l + tanh(a;/2)] (21) 
0±(a;) =±i[l-tanh(a;/2)] (22) 

They are similar to the kink solutions that one finds in the (p^ model, as we have shown in Sec. We can write the 
wall tension in this second model as ^2 = (2/3)is, where tg is the tension value for the BPS solutions of the (j)^ model. 
We illustrate the present model and its topological solutions in Fig. 2. The interesting feature of the second model is 
that its potential is polynomial, of the fourth order type, but it presents two distinct phases, the asymmetric phase 
and another one, symmetric, represented by the vanishing minimum, = 0. 




-1 1 if 



FIG. 1: The potential Vi{4>) (right) given in Eg. 1)10^ . and the respective topological solution (left). For comparison, the dotted 
Une shows the topologial solution of the cj>'^ model. All quantities that appear in these figures are dimensionless. 

The second model seems appropriate to simulate a first order phase transition. However, differently from the 
(j)^ system the present model has the property of having the same squared mass, irrespective of the symmetric or 
asymmetric phases. To see this we use the potential (jllll to obtain 

^ = l-6H + 602 (23) 

so that m2{(j> = 0) = rn2{4' = ±1) = 1- As we shall see below, this fact also shows that U{x) — d?V2/d(lP is symmetric 
respect to x when calculated at any BPS state, so that the investigation to find the first quantum corrections to the 
energy of the BPS solutions follows naturally, cincumventing the intricacy that happens to appear in the standard 
(j)^ model . The form of the potential of Eq. Hll|l shows a symmetry respect to = ±1/2 for in the interval 
< (/> < 1 - see Fig. 2. This feature does not appear in the (j)^ model. 

The classical picture that appears for the second model is somehow similar to the picture of the (p^ model. Thus, 
it may be used as a new model to explore conformational properties of polyethylene (PE), as an alternative to the 
case investigated in Refs. [s^, H^l- Moreover, since its potential contains up to the fourth order power in the field, 
in Sec. we investigate the semiclassical or one-loop corrections, to see how the thermal effects change the classical 
picture. 



C. Model 3 

The third model is given by the potential of Eq. H12() . In terms of the superpotential of Eq. IjlSfl the first order 
equations are 

^^l-|0|-rx^ (24) 
ax 

fx = -'^'^^ 
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FIG. 2: The potential 14(0) (right) given in Eg. llUL and the respective topological solution (left). For comparison, the dotted 
line shows the kink of the 0® model, which is asymmetric. All quantities that appear in these figures are dimensionless. 



There are four degenerate minima for r > 0, at the points vi^2 — (±1,0) and U3.4 = (0, ±l/-\/r). There are five BPS 
sectors, and only one non-BPS sector, which is the sector connecting the minima W3 and W4. The tensions of the BPS 
states are tl = 1 and t| = y/r. 

This model is similar to the model considered in Eq. 0- To find BPS solutions we follow Refs. 0,^3- In the 
sector connecting the minima vi and V2 we have 



tanh(a;/2) 
l + tanh(|x|/2)^ 



Xix) = 



(26) 



The corresponding orbit is a straight line segment joinning the minima vi and V2. We have been unable to find non 
trivial analytical BPS states in this sector, for solutions describing an elliptic orbit, although they appear in the model 
of Eq. ®. 

There are other BPS states, in the sectors that connect the minimum vi or V2 to V3 or V4. There are several cases, 
and we get explicit solutions for the specific value r = 1/2, which form straight line segments connecting vi or V2 to 
V3 or V4. The orbits are described by configurations such that x ~ ±V2(0 T !)• The explicit solutions are, in the 
sectors connecting V2 to or 



-l[l-tanh(|)] 



(27) 



and in the sectors connecting vi to V3 or W4 



X(a;)=±^[l + tanh(|)] 



(x) = i[l + tanh(|)] 



(28) 



(29) 



X(a;)=±^[l-tanh(-: 



(30) 



These solutions are similar to the solutions found in Ref. 52], for the model of Eq. 0. In Fig. 3 we plot one of the 
several BPS states that appear in this case. 

The nontrivial BPS solutions that we found in model 3 can be used to model domain walls solutions that appear in 
models for binary mixtures of Bose-Einstein condensates (BEC) [s^l- As one knows, for different coupling coefficient 
between the two BECs, the mixture exhibit complex spatial structure which may be described by domain walls |54| . 
More importantly, in Ref. [s^ one has realized that these domain walls consist of multicomponent solitons. And these 
multicomponent walls are very much similar to the two-field solutions that we have just found in model 3. A specific 
feature of the BPS state shown in Fig. 3 is that the limit x ^ gives 0(0) > x(0); leading to unequal components in 
the binary mixture. We shall further explore this issue elsewhere. 
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This model is also of interest within the context of of spatially extended systems, where phase fronts separating 
different phase states may appear. The two field model that we are investigating may be used to describe a complex 
order parameter A = (p + ix, A — (p — ix which obeys the Ginzburg-Landau equation 



OA 
'dt 



d\A 



&V 



(31) 



where V = T^(AA) is the potential of Eq. (|12(l . described in terms of the complex order parameter. This equation 
generalizes Eq. to the case of a complex order parameter, and the solutions that we have found in Eqs. (jSEJ, and 
(|?7|l and or and can be used to describe tt, and tt/2 fronts, respectively [39l |. 



^2 ■ 



FIG. 3: A tipical solution that appear from Eqs. I|29|l (thin line) and (dotted line), to illustrate how the fields behave as a 
nontrivial BPS state. All quantities that appear in this figure are dimensionless. 



IV. STABILITY 



The study of linear stability of BPS solutions is of direct interest to high energy physics. The point here is that 
an investigation at the quantum level would require the explicit calculation of eigenvalues and eigenfunctions that 
appear from linear stability. To illustrate the subject, let us investigate linear stability of the BPS solutions obtained 
in the former Sec. IIIII In the case of a single field the investigation requires that one uses 

<j){x,t) = (psix) + ^ rinix) cos{wnt). (32) 

n 

where (f>s{x) stands for the static field, and the remaining terms represent the fluctuations about it. Linear stability 
implies that the fluctuations remains limited in time, so that the set of frequencies {wn} forms a set of real numbers. 
We use Eq. (|32|l into the second order equation of motion (O to obtain the Schrodinger-like equation Hrjn{x) = 
rjn{x), where 

^ = -^ + C^(^) (33) 

and 

U{x) = W^^ - W^W^^^. (34) 

which must be calculated with the static conflguration, the BPS states. 

As we know, we end up with a supcrsymmetric quantum mechanical problem. Thus, we generalize the above 
situation by introducing the pair of Hamiltonians 

H± = —^ + U±ix) (35) 

where 



U±{x)^W^^TW^W^^^. (36) 
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These Hamiltonians can be factorized as 

H± = Ai A± (37) 

where A± are first order operators, given by 

A± = —^±W4,4, (38) 

This factorization ensures that the Hamiltonians H± are positive definite, so that their eigenvalues are all non-negative 
real numbers, thus ensuring linear stability of the classical topologocal solutions. 



A. Linear stability: model 1 

In order to investigate linear stability of the BPS states obtained in the former Sec. IIII AI we use Eqs. (|2Jl, H19|l . 
and H32|l to obtain the Schrodinger-like equation 

+ l-2Six)]rjnix) ^wlvnix) (39) 



In this case the potential is given by 

U{x) = 1 - 26{x) (40) 

Stability of the classical solution H19|l implies that the eigenvalues of the above Schrodinger-like Eq. H39|l should 
be non-negative. But this is indeed the case, because the atractive potential S{x) has only one bound state, and the 
specific form of the potential U{x) = 1 — 26 (x) places this bound state at zero energy. To see this explicitly we follow 
the steps already presented and we factorize the Hamiltonian in the form 



The single bound state is the zero mode, 770(2:), which obeys 



( ^ 




( ^ 




\ dx 




V dx 


x ) 



(41) 



This equation is easily solved to give the normalized zero mode 

= e"'"' (43) 
The quantum mechanical problem is implemented with the partner Hamiltonians. They are 

H± = -^ + U±{x) (44) 

where U±{x) — 1 ± 2 5{x). They are obtained via H± — Al^A±, with the first order differential operators 

d \x\ , , 

dx X 

The bosonic model under investigation can be seen as the bosonic portion of a more general theory, which includes 
fermions via the standard Yukawa coupling. In the model with fcrmions we can search for fermionic zero modes. The 
investigation may follow the lines of Rcf. [Sg ■ 
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B. Linear stability: model 2 

We follow the same steps to investigate stability of the BPS solutions of the second model, given by Eqs. H21|l and 
In this case the Schrodinger-like equation becomes 

("i^ (I)) '^"^^^ " ^IVnix) (46) 



This new Hamiltonian can be fatorized as 



S:-"'-'(f))(-|:-"'°Kf)l <«' 



The zero mode obeys 

d 



tanh(|)j77o(a;)=0 (48) 



dx 

This first order equation can be integrated easily; the normalized zero mode has the explicit form 

1 /3 



'7o(x) = -y-sech^(|) (49) 

We notice that the above Schrodinger-like Eq. H46|l is the equation that appears in the modified Poschl- Teller 
problem, so that the investigation here goes very much like it does in the standard 0^ model - see for instance Ref. 0. 
As a result, besides the above zero mode the present problem has another bound state, at energy wi = V3/2 - recall 
that we are using dimensionless units in the present work. Also, the continuum starts at energy w — 1. 

The quantum mechanical problem is implemented by the partner Hamiltonians H± = Aj_A±, where the first order 
differential operators have the explicit form 



A± = ±- tanh 

dx 



(f) (»' 



C. Linear stability: model 3 

We now direct our attention to the third model. We consider tha classical solution ()26|). In this case we consider 
(j){x,t) as before, and x(x,t) in the form 

X{x,t) = ^ ^n{x) COsiWnt) (51) 
n 

The Schrodinger-like equation splits into two equations, one of them being exactly Eq. (|39|l . and the other is 

d^ 



dx 



^+U{x)]U^)^wlU^) (52) 



where 



_ 8r-(8r + l)sech^(x/2) 
^^"^'-^^ [l-Htanh(N/2)P ^^^^ 

In this case the potential has at least one bound state, the zero mode. From the analytical point of view the problem 
is somehow complicated, and we could not find any explicit solutions. However, we could verify that the number of 
bound states depends on r, and increases for increasing r. As we can see from Eq. H15|) . the parameter r is related 
to the way the two fields interact. Thus, one sees that the number of bound states increases when one increases the 
strength of the interaction between the two fields. This feature also appears in other models, and we shall further 
explore this specific issue in another work. 
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V. THERMAL EFFECTS 



Let us now deal with semiclassical effects. However, instead of calculating the quantum corrections to the energy 
of the classical solutions let us investigate the effective potential in the case of a single real scalar field. There are 
several distinct but equivalent ways of implementing the calculations and we choose to follow Refs. [s^jll^- We get 
the general expression, which is valid at the one-loop level ^3 

(04Mfc'+^"('/')] (54) 

where V"{(j)) = (PV/d^P' , and the metric is now Euclidian. The thermal effects are obtained after changing the one 
loop contribution of Eq. 1)54(1 to 



27rn 



2 

2 



+ V"{(l)) 



(55) 



The teinperature enters the game as /3 = l/T. We make the sum to obtain the finite temperature contribution 

1 f (fik 



The temperature dependent contribution can be used to investigate how the thermal effects enter the game. 
We investigate the second model at finite temperature. We rewrite the potential of Eq. ((TT|l in the form 

V2{4>) = \4>\'m-mf (57) 

Here we have introduced m and A to parametrize the potential, restauring the standard notation that appears with 
Ti = c = \. The parameters m and A are real, and we consider m/A > to allow for spontaneous symmetry breaking. 
We notice that (f) and m has dimension of energy, while A is dimensionless. 

The motivation to study this model is that it has an interesting property, not seem in the (j)^ model. The above 
potential goes up to the fourth order power in the field, and it has minima at the nonzero values |(/)| — to/A and at 
= 0; however, the classical masses corresponding to these minima degenerate to the single (squared) value . 

In the hight temperature limit (T >> m) the one loop effective potential in this case is well approximated by 

K//(0) = \<f'{m - mf + ^ y"(0) (58) 

where V" {(p) = mn? — 6toA|0| + (SX^cfP . We recall that in natural units the temperature has dimension of energy. The 
effective potential has the minima 



l±Jl-—T^\ (59) 





This result allows introducing the critical temperature Tc — to/A, which identifes two distinct behaviors: for T below 
Tc the effective potential supports four minima, at the above values of (pi, z = 1, 2, 3, 4, but for T > Tc there are only 
two minima, at the values (j)± — ±to/2A. The behavior of the effective potential is shown in Fig. 4, where we depict 
two tipical cases, for T = (4/5)Tc and for T = Tc. 

We notice that the high temperature effects are unable to restore the symmetry, which remains broken for T > 
Tc = to/ A, which is greater than to in the weak coupling limit that makes our calculations reliable. Although this 
result goes against conventional wisdow, it has already appeared in other models, as for instance in [sM l6ll | and in 
references therein. 

We investigate the (squared) mass, which can be obtained from the effective potential by the usual procedure. It 
can be written as 

m\T) = V"{^) + ^V""{^)T^ (60) 

Thus, for temperature lower than Tc the masses at the constant configurations pi degenerate to the single value 

to2(T) = TO^ - A^T^ T<Tc (61) 
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For temperature higher than Tc we get 

m'(T) = + yT^, T>T, (62) 

These results show that the effective mass decreases to zero for temperature lower than Tc, and it increases from zero 
for temperature above Tc- The critical temperature identifies the point where the effective mass vanishes, and for 
higher temperature there is no symmetry restoration anymore. 
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FIG. 4: The effective potential Veffi^cf)) as a function of the dimensionless variable = {\/m)(f), for the two tipical values 
T/Tc = 4/5 (right) and T/Tc = 1 (left). The effective potential Veff is dimensionless, given by Veff{<i>) = Veff{4))/X^T^. 

In Fig. 5 we depict the minima of the potential, and the squared mass at finite temperature. There we see that 
(j){T) varies continuosly as the temperature crosses the critical value Tc — m/X, indicating the presence of a second 
order phase transition. 

VI. CONCLUSIONS 

In the present work we have investigated several models described by one and by two real scalar fields. The main 
investigations concern the search for BPS states, that is, for topological solutions that solve first order differential 
equations. These solutions minimize the energy to the Bogomol'nyi bound, which is given solely in terms of the 
superpotential, and the asymptotic value of the corresponding field configurations. 

The search for topological solutions is done at the classical level, and we have payed special attention to models 
introduced in Sec. ^1 some of them further explored in Sec. IIIII These investigations have shown that model 1 and 
model 2, defined by the potentials of Eqs. H1U|) and Ijlll) . respectively, are classically similar to the standard field 
theory models known as the 0^ and the cffi models. The other model is a two-field model, defined by the potential of 
Eq. H12|) . It is similar to a model first investigated in Ref. |4§l] . 

The classical investigations that we have developed are of interest in applications to nonlinear science, as for 
instance in the line of the work presented in Refs. |3^ |3^ where one uses field theory models to mimic 
nonlinear interactions in polimeric chains. For instance, the model 1 may be used as an alternative to the model 
that is standardly considered to mimic the polyacethylene (PA) chain, where Peierls instability appears due to single 
and double bound alternation in carbon atoms along the chain. The standard scenario leads to the very nice picture 
in which the distance carbon-carbon is maped to the 0^ model with spontaneous symmetry breaking. The need for 
spontaneous symmetry breaking is to reproduce the two degenerate states, which describe single-double and double- 
single alternations in the trans or zig-zag PA chain. Of course, this picture can be more interesting if one adds 
fermions to the system, via the standard Yukawa coupling. The alternative that we propose is to mimic the PA chain 
with model 1, which is very much similar to the 0"* model at the classical level. We hope to explore this and similar 
ideas in the near future. 

In order to further explore the analogy between standard field theory models and the models here introduced, we 
have also calculated the thermal effects at the one-loop level. In particular, we have explored the second model, 
working out the one-loop finite temperature correction to see how the semiclassical contributions may change the 
classical picture of the model. As we have shown, the finite temperature effects add to the classical potential of model 
2 and change it in a way such that it allows for a second order phase transition to take place at high temperature, 
although the symmetry is never restored by the thermal corrections. 
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FIG. 5: The (dimensionless) minima 4>{T) [(f> = (A/m)0] of the effective potential (right), and the effective mass vn?{T) (left), 
depicted as a function of the temperature. Here we are using natural units, thus the mass m and the temperature T have 
dimension of energy, and are measured in electron- Volt. The dotted-line regions in both figures are depicted to inform that our 
results are obtained within the high temperature approximation, and they may not be valid in low temperatures. 



The potential of the second model is polynomial, and goes up to the fourth order power in the field. In this sense 
it reminds us of the 0"* model. However, it admits two degenerate but different phases, the symmetric phase that is 
governed by the vanishing minimum, and the asymmetric phase which is described by the two non- vanishing minima. 
In this sense it is similar to the cff' model. However, at the semiclassical level the thermal corrections change the 
behavior of the model in a way such that it is very specific, reminding us of neither the (j)^ nor the 0^ models anymore. 

A direct motivation that follows from the present investigations concerns the inclusion of fermions, to see how 
the fermions change the scenario we have just obtained. Another motivation concerns generalization of the field 
theory models to the case of complex fields, which give rise to models engendering the continuum U{\) symmetry. 
In this new scenario the models admit the presence of vortices, global and local, depending on the gauging of the 
global symmetry that appears when one changes the real field to a complex one. Furthermore, the two-field model 
investigated in Sec. IIIII is of direct use to describe interactions between two EEC's, and the corresponding vector 
solution describes the interface between the interacting condensates, as a good alternative to the recently proposed 
model of Ref. 55| since in our model the center of the interface is asymmetric, containing different quantities of each 
one of the two condensates that form the interface. This is more realist than the case studied in where the center 
of the interface contains equal portions of each one of the two condensates. All the models here introduced are also 
of interest to pattern formation in spatially extended, periodically forced systems, governed by the Ginzburg-Landau 
equation. This is a new route, different from the one proposed in [sj], and further investigated in |32ll33l l33|. which 
deals with solitons in ferroelectric cristals, in polyethylene, and in Langmuir films. The Ginzburg-Landau equation 
is appropriate for investigating fronts and interfaces, and their contributions to pattern formation. These and other 
specific issues are presently under investigation, and we hope to report on them in the near future. 

We thank M. Bahiana, P.L. Christiano, J.R.S. Nascimento, R.F. Ribeiro, and C. Wotzasek for interesting discus- 
sions, and CAPES, CNPq, PROCAD and PRONEX for financial support. 
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